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Abstract 

Given a finite sequence of times < t\ < ■ ■ ■ < , we construct an example of 
a smooth solution of the free nonstationnary Navier-Stokes equations in R d , d = 2,3, 
such that: (i) The velocity field u(x,t) is spatially poorly localized at the beginning of 
the evolution but tends to concentrate until, as the time t approaches t\ , it becomes 
well-localized, (ii) Then u spreads out again after t\ , and such concentration-diffusion 
phenomena are later reproduced near the instants £2 , £3 , 

1 Introduction 

One of the most important questions in mathematical Fluid Mechanics, which is still far from 
being understood, is to know whether a finite energy, and initially smooth, nonstationnary 
Navier-Stokes flow will always remain regular during its evolution, or can become turbulent 
in finite time. 

As a first step toward the understanding of possible blow-up mechanisms, it is interesting 
to exhibit examples of smooth and decaying initial data such that, even if the corresponding 
solutions remain regular for all time, "something strange" happens around a given point {xq, to) 
in space-time. This is the goal of the present paper. 

Our main result is the construction of a class of (smooth) solutions to the incompressible 
Navier-Stokes equations such that, in the absence of any external forces, the motion of the 
fluid particles tends to be more concentrated around x$ , as the time t approaches to • This 
corresponds precisely to the qualitative behavior that one would expect in the presence of 
a singularity, even though such "concentration of the motion" is not strong enough to imply 
their formation. 

1.1 Statement of the main result 

For a fluid filling the whole space M. d , d > 2 , the Navier-Stokes system can be written as 

' d t u + PV-(u®m)=Au 
< V -u = 
k u(x, 0) = a(x), 

* Keywords: Navier-Stokes, asymptotic profiles, asymptotic behavior, far-field, flow map, analyticity, sym- 
metry, spatial spreading, lower bound, pointwise estimates, decay rate, singularities. 
2000 Mathematics subject classification: 76D05, 35Q30 

^Universite de Lyon, Universite Lyon 1, CNRS UMR 5208 Institut Camille Jordan, 43, bd. du 11 novembre 
1918, F - 69622 Villeurbanne Cedex, France. E-mail: brandolese@math.univ-lyonl.fr 



where u = (u\, . . . , Ud) , a is a divergence-free vector field in M. d and P = Id — VA _1 div is 
the Leray-Hopf projector. 

Because of the presence of the non-local operator P a velocity field that is spatially well 
localized (say, rapidly decaying as \x\ — > oo) at the beginning of the evolution, in general, 
will immediately spread out. A sharp description of this phenomenon is provided by the two 
estimates (2) below. In order to rule out the case of somewhat pathological flows (such as two 
dimensional flows with radial vorticity, or the three-dimensional flows described in [3], which 
behave quite differently as \x\ — > oo if compared with generic solutions), we will restrict our 
attention to data satisfying the following mild non-symmetry assumption (for j, k = 1, . . . , d): 

3j*k: /(a^M^O, or J \ j(x f <l X * j <fa. (1) 

Then, for sufficiently fast decaying data, we have, for \x\ > ^= (see [4]): 

r?i(i)W~ d_1 < KM)| < m{t)\x\- d ~\ |f| G S d - 1 \S a , (2) 

these estimates being valid during a small time interval t G (0,ti). Here C, t\ > and 771 
and rj2 are positive functions, independent on x, behaving like ~ Cj t as t — > (j = 1,2). 
Moreover, S^ 1 denotes the unit sphere in M. d and the subset S a of S rf_1 represents the 
directions along which the lower bound may fail to hold: the result of [4] tells us that E a can 
be taken of arbitrarily small surface measure on the sphere. In other words, the lower bound 
holds true in quasi-all directions, whereas the upper bound is valid along all directions. 

Moreover, the upper bound will hold during the whole lifetime of the strong solution u 
(see [12], [13]), whereas the lower bound is valid, a priori, only during a very short time 
interval. The main reason for this is that the matrix (f UjUk{x,t)dx) is non-invariant during 
the Navier-Stokes evolution, in a such way that even if the datum satisfies (1), it cannot be 
excluded that at later times the solution features some kind of creation of symmetry, yielding 
to a better spatial localization and, after t > t\ , to an improved decay as \x\ ^00. (We refer 
e.g. to [2], [3], [7], [8], for the connection between the symmetry and the decay of solutions). 

The purpose of this paper is to show that this indeed can happen. We construct an 
example of a solution of the Navier-Stokes equations, with datum a G 5(IR rf ) (the Schwartz 
class), d = 2, 3, such that the lower bound 

\u(x,t)\ >Vi(t)\x\- d -\ (3) 

holds in some interval (0, t±) , but then brakes down at t\ , where a stronger upper bound can 
be established. This means that the motion of the fluid concentrates around the origin at 
such instant. Then the lower bound (3) will hold true again after t\ , until it will break down 
once more at a time t<i > t± . This diffusion-concentration effect can be repeated an arbitrarily 
large number of times. 

More precisely, we will prove the following theorem. 

Theorem 1 Let d = 2, 3, let < t\ < • • • < tjy be a Unite sequence and e > 0. Then there 
exist a divergence-free vector field a G and two sequences (t[, . . . , t' N ) and (t*, . . . , t* N ) 

such that the corresponding unique strong solution u(x, f) of the Navier-Stokes system sat- 
isfies, for all i = 1, . . . , N and all \x\ large enough, the pointwise the lower hound 
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and the stronger upper bound 

\u(x,t*)\ <C\x\- d ~ 2 , 

for a constant C > independent on x and a constant c w independent on \x\, but possibly 
dependent on the projection oj = ^ of x on the sphere, and such that c w > for a.e. 

lo G S^" 1 . Moreover, t\ and t* can be taken arbitrarily close to U : 

— ti\ < e and \t* — tj| < e, for i = 1, . . . ,N. 

Remark 1 The initial datum can be chosen of the form a = curl(^) , where ip is a linear 
combination of dilated and modulated of a single function (or vector field, if d = 3) (f> € <S(IR d ) , 
with compactly supported Fourier transform. 

Roughly speaking, our construction works as follows: we look for an initial datum of the 
form a = curl(^) , where 

d(N+l) 

ip(x) = ^2 \j 5 d / 2 <p(5x)cos(oij -x). 

3=1 

The unknown vector a = (cci, . . . , awjv+i)) G M d2 ( JV+1 ) of all the phases ctj G lR d will be 
assumed to belong to a suitable subspace V C 

R d 2 (jv+i) of dimension ^(N + 1) in order 
to ensure, a priori, some nice geometrical properties of the flow. Such geometric properties 
consist of a kind of rotational symmetry, similar to that considered in [2], but less stringent. 
In this way, the problem can be reduced to the study of the zeros of the real function 

J (u\U2)(x, s) dx ds. 

By an analyticity argument, this in turn is reduced to the study of the sign of the function 

1 1— ► J J (e sA aie sA a2)(x, s) dx ds. 

This last problem is finally reduced to a linear system that can be solved with elementary 
linear algebra. 

The spatial decay at infinity of the velocity field is known to be closely related to special 
algebraic relations in terms of the moments J x°curl(u)(x, t) dx of the vorticity curl(tt) of the 
flow, see [6]. Thus, one could restate the theorem in an equivalent way in terms of identities 
between such moments for different values of a £ N^, which are satisfied at the time t* but 
brake down when t = t\ . 

1.2 A concentration effect of a different nature 

The concentration-diffusion effects described in Theorem 1 genuinely depend on the very 
special structure of the nonlinearity PV ■ (u®u), more than to the presence of the Am term. 
Even though this result is not known for inviscid flows yet, it can be expected that a similar 
property should be observed also for the Euler equation. 

On the other hand, the Laplace operator, commonly associated with diffusion effects, can 
be responsible also of concentration phenomena, of a different nature. For example, it can 
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happen that a(x) is a non decaying (or very slowly decaying) vector field, but such that 
the unique strong solution u(x, t) of the Navier-Stokes system have a quite fast pointwise 
decay as \x\ — > oo (say, ~ This is typically the case when a has rapidly increasing 

oscillations in the far field. We will discuss this issue in Section 3. Though elementary, the 
examples of flows presented in that section have some interest, being closely related to a 
problem addressed by Kato about strong solutions in L p (M. d ) when p < d, in his well known 
paper [9]. 



1.3 Notations 

Troughout the paper, if u = (ui, . . . ,Ud) is a vector field with components in a linear space 
X , we will write u G X , instead of u G X d . We will adopt a similar convention for the 
tensors of the form u®u. We denote with e iA the heat semigroup. 

Let B(0, 1) be the unit ball in R d and <p G 5(M d ) a function satisfying 



4>eC^(R d ), supp^C 5(0,1), radial, > 0, J\4>\ 2 = l/d. (4) 
Our definition for the Fourier transform is (/>(£) = J (f>(x)e~ l ^' x dx . Then we set 

$ S (0 = * > 0. (5) 



(6) 



Next we define the orthogonal transformation : M. d — ► M. d , by 

a = (cc2,oti), if a = (a±, 02) G M 2 , 

5 = (062, as, ai), if a = («i, «2, 03) G M 3 . 

We define the curl(-) operator by 

curl^ = (-d 2 ,di)ip, if^:K 2 ^M 

and by 

($2^3 ~ ^3^2 \ 
^3^1-91^3 if^:IR 3 ^lR 3 - 

The notation f(x,t) = Ot(\x\~ a ) as \x\ — > 00 means that / satisfies, for large \x\ , a bound 
of the form |/(x,i)| < ^(t)^) - ", for some function A locally bounded in M + . 

We shall make use of the usual Kronecker symbol, Sj t k = 1 or , if j = k or j / k . 



2 Nonlinear concentration-diffusion effects 

2.1 The analyticity of the flow map 

In this subsection we recall a few well known facts. 

Let B be the Navier-Stokes bilinear operator, defined by 

B(u, v)(t) = - [ e {t -^ A FV • (u ® v)(s) ds. 
Jo 
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Then the Navier-Stokes equations can be written in the following integral form 

u = uq + B(u, u), uo = e tA a, div(o) = 0. (7) 

Even though in the sequel we will only deal with "concrete" functional spaces, the prob- 
lematic is better understood in an abstract setting: we will present it as formulated in the 
paper by P. Auscher and Ph. Tchamitchian [1]. Let T be a Banach space, uq £ T and let 
B : T x T — > T be a continuous bilinear operator, with operator norm . 

Let us introduce the nonlinear operators T k : T — > T ', k = 1,2..., defined through the 
formulae 

Tl = Idjr 

fc-1 

T k (v) = J2 B ( T l( v ^ T k-l( v »' k ^ 2 ' 
l=i 

Then the following estimate holds: 

\\T k (u )y<^k- 3 / 2 (4\\B\\ \\ Uo y) k . 

Moreover, T k is the restriction to the diagonal of T k = T x • • • xJ 7 of a k -multilinear operator 
: T k -> T . 

Under the smallness assumption 

||«o||^ < 1/(4||B||), (8) 

the series 

oo 

tf(u ) = Y, T k(u ), (9) 

k=l 

is absolutely convergent in T and its sum \l/(uo) is a solution of the equation u = uq+B(u, u) . 
Furthermore, ^(uq) is the only solution in the closed ball -Bjf(C), 2pjj[)- The proof of these 
claims would be a straightforward application of the contraction mapping theorem under the 
more restrictive condition ||«o||.F < V(4||-B||). The proof of the more subtle version stated 
here can be found in [1], [11]. 

Coming back to Navier-Stokes, in the proof of Theorem 1 we will need to write the solution 
of the Navier-Stokes system as u = «£(a), where 

oo 

d>(a)(t)^^T fc (e' A a), (10) 
k=i 

the series being absolutely convergent in C([0, oo), L 2 (M. d )) . There are several ways to achieve 
this, and one of the simplest (which goes through in all dimension d > 2) is the following: we 
choose T as the space of all functions / in C([0, oo), L 2 (U. d )) , such that < oo, where 

\\fy = e SS sup Xtt (l + \x\) d+1 \f(x,t)\ + esssup^a + ^+^l/OM)!. (11) 
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The bicontinuity of the bilinear operator B is easily proved in this space T (see [12]). Indeed, 
one can prove this only using the well known scaling relations and pointwise estimates on the 
kernel F(x,t) of the operator e* A Pdiv: 

F(x, t) = t-( d+1 V 2 F{x/Vi, 1), \F(x, 1)| < C(l + \x\)- d -\ 

We can conclude that there is a constant % > 0, depending only on d, such that if 

||e' A a||^ < Vd (12) 

then there is a solution u = <J?(a) £ J 7 of the Navier-Stokes equations such that the se- 
ries (10) is absolutely convergent in the ^"-norm. The absolute convergence of such series 
in C([0, oo), L 2 (R. d )) is then straightforward under the smallness assumption (12). 

The finiteness of ||e tA a||jF can be ensured, e.g. by the two conditions f |a(x)|(l + \x\) dx < 
oo and ess sup^j^l + |x|) d+1 |a(x)| < oo. The smallness condition (12) could be be slightly 
relaxed, see [4]. 



2.2 The construction of the initial datum 

This section devoted to a constructive proof of the following Lemma. 

Lemma 1 Let d = 2,3, and e > 0. Let also N G N and < t\ < • • • < tjy be a finite 
sequence. Then there exists a divergence-free vector field a = (a±, . . . , ad) £ 5(M d ), such that 

a(x)=a(x), xeR d (13) 

(see Section 1.3 for the notations) and such that the function E(a)(t) : M + — > R, defined by 

E(a){t) = ~ J^J e sA a 1 (x)e sA a 2 (x)dxds, (14) 

changes sign inside (tj — e, tj + e) , for i = 1, . . . , N . 

Proof. It is convenient to separate the two and three-dimensional cases 
The case d = 2 . We start setting, for each a£l 2 , 

MO = HZ - a) + + a) - 0(£ - 5) - ^ + 5), (15) 

for some (j) G 5(]R 2 ) satisfying conditions (4). Next we introduce the divergence-free vector 
field a a (x) , through the relation 

MO^tf). (16) 



^lV'a(C) 

Note that a Q G 5(M 2 ) is real- valued (because ty a is real- valued and such that Y> a (£) = 
V'aC - 0) an( l satisfies the fundamental symmetry condition 

a a (x) = a a (x). (17) 

Next we define tp s , as before, by simply replacing (f> with (f> s in the corresponding defini- 
tions. The vector fields a 5 a will be our "building blocks" of our initial datum. 
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Applying the Plancherel theorem to the right-hand side of (14) and using the symmetry 
relations ^ Q (£) = $ a (Q\ = |V>a(f)|, we get 



E(a a )(t) = j(l-e-^)^\MOM 



(18) 



= 2/ (i-e-^)^|v,,iir</i. 
Jti>\to\ 



Prom now on, the components of a G I 2 will be assumed to satisfy the following conditions 

h > M (19) 

This guarantees that for a sufficiently small 5 > (i.e. when ai > |ck2 | + we have 

E{ai){t) = 2j(l-e-W\ 2 )^\fi{t-a)\ 2 dt. (20) 

If we set 

^ (t) = (l_ c -*M a )2^ (21) 

then we immediately obtain 

E(a s a )(t)^E%*(t), as<5^0 

uniformly with respect to t > 0. 

We now associate to (t±, . . . , tjy) two more sequences (to be chosen later) (Ai, . . . , Atv + i) G 
and (a±, . . . , ayv+i) G M 2 ( Af+1 ) , where aj = (aj,i, 0^,2) • First we require that the 
components a^i and 0^2 satisfy condition (19) for all j = 1, . . . , N + 1 and that aj 7^ a^/ , 
for j / j' and j,f = 1, . . . , N + 1 . This second requirement ensures that the supports of a* 
and a^., are disjoint when 5 is sufficiently small. 

We now consider the initial data of the form 

N+l 

a s (x)=J2 X AW- ( 22 ) 



Owing to the condition on the supports of a* , we see that for 5 > small enough 



N+l 

E(a s )(t) = Y,^E(a s aj )(t). (23) 
j'=i 



Let 

\2 



A.- CKo 1 flj 2 nl 12 

= 3 , l9 and A,- = e- 2 l^l . (24) 

Thus, as (5 -»■ 0, we get E(a s )(t) -► £ a PP(t), uniformly in [0,oo), where 

JV+l 

2^(t) = 5>,-(l-A$). (25) 
j'=i 
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Let us observe that 



d£ app 
dt 



N+l 
3=1 



(26) 



We want to determine (Ai, . . . , \n+i) and (cci, . . . , ayv+i) in a such way that E app (t) vanishes 
at ti,...,t/v, and changing sign at those points. This leads us to study the system of iV 
equalities and A/" 'non-equalities', 



' E app (U) = 
d£ app 



dt 



■(ti) 7^ 0. 



i = l,...,N. 



(27) 



Let us choose \oij\ 2 



7J, for an arbitrary 7 > (this choice is not essential, but will 
ilify the calculations) and set Tj = e~ 2lU . Recalling (24), we get A*? = T{ and 
2jj . In order to study the system (27), we introduce the (N + l) 2 -matrix 

(\-T x 1-7? ••• l-7f +1 \ 



M :-- 



1-T N l-Tl 



N 



1 - 



(28) 



V Ti 27? ••• (N + l)Tf r+1 J 
We claim that detM / 0. Indeed, by an explicit computation, 



N N 

det(M) = -T l (l-T 1 )Yl(l-T l )H(T 1 -T i ) ]J {T v -T t ). 

l<i<i'<N 



i=l 



i=2 



Recalling that Tj = e~ 2jti G (0, 1) and that ti 7^ tj/ proves our claim. The above formula can 
be checked by induction. Otherwise, one can reduce M after elementary factorizations to a 
Vandermonde-type matrix (see [10] for explicit formulae on determinants). 
Then, for any c^O, the linear system with unknown fj, = (fii, . . . , fJ,N+i) , 



AA 



w 



(29) 



has a unique solution /x* 6 M Ar+1 , /x* ^ 0. By our construction, the function .E app (t) ob- 
tained taking /1 = /x* satisfies the iV-equations and the first 'non- equality' of the system (27). 
More precisely, we get (d/dt)£' app (ti) = 0/(27) / 0. The other N — 1 'non-equalities' of the 
system (27) are then automatically fulfilled. Indeed if, otherwise, we had (d/dt)£' app (tj) = 0, 
for some i = 2, . . . , N , then the matrix obtained replacing in (28) the last line with 

Ti 27? ■■■ (N + l)T t N+1 

would have been of determinant zero, thus contredicting our preceding formula for det(M). 

By conditions (24), for all j = 1, . . . , N + 1 , the real number /x* defines (in a non- 
unique way) a real A* and a vector (a* 1; a* 2 ) with components satisfying (19), such that 
;A. A.v.;i / (0 ()). 
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We now consider the initial data obtained from formula (22), choosing Xj = X* and 
aj = a* for j = 1, . . . , N + 1 . With this choice, the corresponding function E app (t) changes 
sign in each neighborhood of tj . 

By the uniform convergence of E(a%)(t) to E app (t) as 5 — ► 0, we see that if 5 > is small 
enough then E(a°) changes sign in the interval (tj — e, U + e) , for i = 1, . . . , N . 

The conclusion of Lemma 1 in the two dimensional case now follows. 

The case d = 3 . We only indicate the modifications that have to be done to the above 
arguments. Let us set, for a£l 3 , 



Oa(0 

where ip is defined by 



(«6-^iM)l ( :3 °) 



^a(o =0(e-«) + 0(e-5) + ^-5) 

+ 0(£ + a) + 0(£ + 5) + + 5). 

In this way, a G 5(]R 3 ) is a real valued divergence-free vector field with the rotational sym- 
metry (13). As before, we can define also the rescaled vector field a s , for any 5 > 0. 
The conditions to be imposed on the components of a = (cci, ct2, as) are now 

min(Q 2 ,a3) > £*i~, «2 / ®3, (31) 

where = max(ai,0). Geometrically, the inequality in (31) corresponds to cutting M 3 into 
six congruent regions, that can be obtained from each other through the orthogonal transforms 
a«a and a 1— ► —a, and then selecting one of these regions. If 5 > is small enough, in a 
such way that a + B(0,5) is contained in the region T C M 3 defined by (31), then we get, 
recalling the definition (14) of E(a a ) , the three-dimensional counterpart of (20): 

E{a a )(t) = 3/(1 - e-^ 2 ) ^-^-^ \?(C ~ «)| 2 di. 

(Here one first applies Plancherel Theorem, then the integral over £ 6 M 3 is rewritten, because 
of the symmetries of a a , as 6 J r . . . c/£). We have, for all i > and as <5 — ► 0, 



£7(a*)(t) - £r (t) = (1 - e- 2 *H 



2 (ai - a3)(«2 - "3) 



lap 



For aj = (ct;,i, <Xj,2, a j,3) G M 3 , ,7 = 1, . . . , JV + 1 , with components satisfying (31), we now 
set 

A 2 (aj-,1 - a ii3 )(aj,2 - ttj-,3) 
^ = ^ j^-p . (32) 

As in the two-dimensional case, it is possible to choose the phases aj in a such way that 
\ a j\ 2 = 13 i where 7 > is arbitrary. This leaves unchanged the definitions of Aj and Tj. 

Now solving, as before, the linear system (29), shows that it is possible to construct an 
initial datum of the form 

jV+l 

3=1 
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such that E(a)(t) has a non-constant sign inside the intervals (U — e, tj + e), for i = 1, . . . , N . 
Lemma 1 is now established. 

□ 

2.3 End of the proof of Theorem 1 

We are now in position to deduce from Lemma 1 and the facts recalled in Section 2.1 the 
conclusion of Theorem 1. 

Step I. Constructing a solution such that t f*J uiu 2 (x, s) dx ds changes sign near t\, 

t 2 , • • • , <JV ■ 

Let us consider the initial datum a constructed in Lemma 1. If necessary, we modify a by 
multipliying it by a small constant rjo > in order to ensure that the corresponding solution u 
of the Navier-Stokes system is defined globally in time. Without loss of generality, we can 
and do assume 770 = 1 • 

With the same notation of Section 2.1, let K{a) : R+ -> R dxd , 

K(a)(t) = j J ($(a) <g> $(a)) (x, s) dx ds. (33) 

Let T be the space defined in Section 2.1, and normed by (11). For 77 > small enough we 
can write, by expansion (10), 

00 

K(r,a)(t) = J2v k S k (e tA a)(t) (34) 

k=2 

where S k ■ J 7 —> R dxd is the restriction to the diagonal of T k = T x • • • x T of a A; -multilinear 
operator defined from T k to M dxd . 

The embedding of T C C([0, oo),L 2 (R d )) implies that, for all t > 0, the series in (34) is 
absolutely convergent. But, 

S 2 (e tA a)(t) = J^J {e sA a®e sA a){x,s)dxds, 

so that, according to the notations of Lemma 1, 

(S 2 (a)) h2 (t) = -E(a)(t). 

For any time t such that E(a)(t) > (respectively, E(a)(t) < 0), we can find a small rj t > 
such that the solution $(77 a) of the Navier-Stokes system starting from 77a satisfies, for all 
< 77 < rjt , 

K(rja)i j2 (t) < 0, (respectively, K(rja)i^(t) > 0). 

By Lemma 1, this observation can be applied for t belonging to two suitable finite sequences 
of times (that can be taken arbitrarily close to (t\, . . . , ijv) in the M^-norm). We conclude 
that if 77 > is small enough then K(rj a)\^{t) changes sign in (ij — e, U + e) , for i = 1, . . . ,iV. 
In particular, because of the continuity of the map t 1— > K(rja)(t), there is a point t* inside 
each one of these intervals where K(r] a)i^(t*) = 0. 

Step 2. Analysis of the orthogonality relations JqJ (ujU k )(x, s) dx ds = c(t)6j jk . 
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Let u = $(77 a) be the solution constructed in the first step. Notice that u(x,t) = u(x,t) for 
all t > 0, i.e., the condition a(x) = aix) propagates during the evolution. Indeed, this is 
simple consequence of the invariance of the Navier-Stokes equations under the tranformations 
of the orthogonal group 0(d) and the uniqueness of strong solutions. Thus, fufdx = 
f u?,dx. When d = 3, such integrals equal, of course, f u^dx, and we have also J u\U2 dx = 
J U2U3 dx = J U3U1 dx . We deduce that, for all t > (we denote here by / the dx d identity 
matrix) : 

3c(t): K( V a)(t) = c(t)I if and only if K( V a) h2 (t) = 0. (35) 

(This equivalence is no longer valid for d > 4. For a proof of the theorem in the higher 
dimensional case one should consider flows invariant under larger discrete subgroups of 0(d)). 

Step 3. The far-field asymptotics of the velocity field. 

By the result of [4] (see Theorem 1.2 or Theorem 1.7, applied in the particular case of the 
datum rja G 5(M rf )) we know that, for all t > 0, 

u(x,t) = V x U(x,t) + Ot(\x\- d ~ 2 ), as|x|^oo, (36) 

where, 

n (^) = ^EGIf - W*) ■ (37) 

h,k V 1 1 11 7 

and 7^ 7^ is a constant. 

Concerning the first term V X H on the right-hand side of (36), for each fixed t > two 
situations can occur. Either the function x 1— > V^II^, t) is identically zero, or this function 
is homogeneous of degree exactly — d — 1. But, for all fixed t > (see [4], Proposition 1.6), 

VJI(-,t) = if and only if 3 c(t) G R: K(r)a) h , k (t) = c(t)6 h , k . (38) 

Combining conditions (35)-(38) with the asymptotic profile (36), we deduce from the 
analysis we made in Step 1, the upper bound 

\u(x,t*)\ < C\x\~ d ~ 2 , i = l,...,N 

and the lower bound, 

> cjxl" 6 *" 1 , i = l,...,iV, j = l,...,d, 

for all |x| large enough and some points distant less than e from t j . Here, C > is 
independent on x and is independent on |x| , but will depend on the projection lo = ^ 

of x on the sphere In fact, we can take c w > 0, unless <? x .II(cj,£^) has a zero at the 

point lo G S d_1 . But one deduces from (37) that the zeros of d Xj U.{u},ty are exactly the 
zeros on the unit sphere of a homogeneous polynomial of degree three. Therefore, > for 
almost every lo G § rf_1 . 

Theorem 1 is now established. 

□ 
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2.4 Explicit examples 

Let us exhibit an explicit example of our construction, in the simplest case N = 1. We set, 
for d = 2 , and a function <j> satisfying conditions (4) , 



a K x ) = V ( di 



lS (x)[ v / 3cos(v / 3~£i + ^2) — v / 3cos(a;i + \/3x2) 

+ y/2 cos(\/6xi — V2X2) — V2 cos(— \/2a:i + \/6:E2)) 



This expression is obtained taking, in Lemma 1, Ti = ^ , ai = (\/3, 1) , 02 = (\/6, — \/2) and 
7 = 4, and observing that Af = |A|. If and 5 are positive and small enough, then the 
solution u(x,t) of Navier-Stokes starting from a concentrates/diffuses when t ~ |log(2). 

Three-dimensional examples are obtained in a very similar way: choose again T\ = ^ , 
next a\ = (0, l,\/3), «2 = (0, V6, V2) ■ The relation between the coefficients Ai and A2 is 
now A 2 = -^A 2 . This leads us to introduce the function 



f\ X ) EE Cf> 5 (x) 



^cos(a;2 + VSx^) + cos(xi + V3X2) + cos(\/3~£i + 2:3)^ 
+ v / 2^cos(v / 6a;2 + V^x^) + cos(\/6:ri + V2X2) + cos(\/2:ei + \/6x3)^ 



(39a) 

where (5 > and satisfies conditions (4). If a is the vector field 

a = r ] curl(f s ,f 5 ,f), (39b) 

with \r]\,5 > small enough, then the solution arising from a concentrates/diffuses, as before, 
around t ~ | log(2) . 

Remark 2 The smallness of 77 was important for the applicability of our analyticity argu- 
ment. However, it would be interesting to know whether three-dimensional data as those 
constructed in (39a)-(39b), with large coefficients 77, still feature some kind of concentration 
effects in finite time. 



3 Linear concentration-diffusion effects 

In this section we would like to give an elementary example of a Navier-Stokes flow featuring 
a concentration effect of a quite different nature. A similar effect can be observed if we replace 
the Navier-Stokes nonlinearity by a more general term (not necessarily quadratic). Let us 
focus, for example, on the three-dimensional case. 

Consider the two initial data in R 3 (here rj 7^ is a constant) 

a(x) =r/(-a 2 [log(e+|x| 2 )- 1 ] , 8 1 [log(e + M 2 )" 1 ] , 0) (40) 

and 

a(x) = a(x) sin(|x| 2 ). (41) 
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In both direct computation shows that such data are divergence-free and belong to 

L 3 (M 3 ) , but not to LP(R 3 ) for 1 < p < 3. 

In his well-known paper [9], Kato could obtain, for a solution u of the Navier-Stokes 
system emanating from a, local and global results implying u(t) £ L q (M. 3 ) for t > 0, with 
q > 3, under the assumption a £ L 3 (IR 3 ). He also obtained results implying u(t) £ L p n 
L 9 (R 3 ) , always for <? > 3, under the stronger assumption a £ L p n L 3 (M 3 ) , with 1 < p < 3. 
Immediately after stating his theorems, Kato observed the following 

Remark 3 (see [9]) "The spatial decay expressed by the property u(t) £ L q (W 3 ) is of 
interest. Note that q is restricted to q > 3. W^e were aft/e to giwe no results for q < 3 under 
the [only] assumption a £ L 3 (R 3 )". 

The solutions arising from the data d and a (both uniquely defined in C([0, T), L 3 (M 3 )) , 
for some < T < oo a priori depending on 77) show that when one assumes only a 6 L 3 (IR 3 ) 
then "everything can happen". Indeed, it is not difficult to prove the following claims: 

1. The solution u(t) arising from the datum a does not belong to L 9 (IR 3 ) , for any t £ (0, T) 
and any q £ [1, 3) . 

2. The solution u(t) arising from the datum a does belong to L 9 (]R 3 ) for all t £ (0,T) and 
all q £ [1,3). 

Therefore, in the latter case, the solution enjoys some kind of spatial concentration effect. 
However this effect does not rely on special geometric features of the flow, but only on the 
oscillatory character of a. 

The first claim is immediate, because whenever \a(x)\ < C(l + H) -1 and |Va(x)| < 
C(l + |x|) -2 , then arguing as in [5, Sec. 4] one obtains for t > an estimate of the form 
\e tA a(x) — a(x)\ < C(t)(l + \x\)~ 2 . Then, from the equation u = e tA a + B(u, u) (the notation 
is as in Section 2.1), the estimate \u(x,t) — a(x)\ < C(t)(l + \x\)~ 2 . Applying this observation 
to d yields the conclusion. 

Let us sketch the proof of our second claim. The faster and faster oscillations of a imply 
that, for t > and \x\ large enough, 

|e' A a(x)| < C m (t)|xr m_1 log- 2 |x|, m = 0,1,2... 

where the functions C m (t) are locally bounded in (0, 00) and satisfy C m (t) ~ t" m/2 as t -> 0. 
In particular, ||e iA a|| 3 /2 < C\(t) . Let us iterate the integral Navier-Stokes equation: 

u (t) = e tA a + B(e tA a, e tA a) + 2B{e tA a, B(u, u)) + B(B(u, u),B(u, u)). (42) 

Now applying elementary Holder and Young inequalities in the right-hand side of (42) (we can 
freely use here that ||u(£)||3 is bounded in (0,T)), we get, for all t £ (0,T), ||it(i)||i < C(t) 
with C(t) ~ t~ l as t ->■ 0. By interpolation one obtains that u{t) £ L 9 (R 3 ) in (0,T) for 
1 < P < 3, with an estimate on the blow-up of the L^-norm as t — > 0. 

Let us observe that for < t < T, even though u(t) decays pointwise much faster than a, 
as |x| ^00, the decay of u(t) does not exceed that of |x| -4 , accordingly with the limitations 
on the spatial localization described in the introduction. 
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